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Abstract 

We study the phase transition of asymmetric hadronic matter to a quark- 
gluon plasma within the framework of a simple two-phase model. The analysis 
is performed in a system with two conserved charges (baryon number and 
isospin) using the stability conditions on the free energy, the conservation laws 
and Gibbs’ criteria for phase equilibrium. The EOS is obtained in a separate 
description for the hadronic phase and for the quark-gluon plasma. For the 
hadrons, a relativistic mean-field model calibrated to the properties of nuclear 
matter is used, and a bag-model type EOS is used for the quarks and gluons. 

The model is applied to the deconfinement phase transition that may occur in 
matter created in ultra-relativistic collisions of heavy ions. Based on the two- 
dimensional coexistence surface (binodal), various phase separation scenarios 
and the Maxwell construction through the mixed phase are discussed. In 
the framework of the two-phase model the phase transition in asymmetric 
matter is continuous (second-order by Ehrenfest’s definition) in contrast to 
the discontinuous (first-order) transition of symmetric systems. 


*Present address: TRIUMF, 4004 Wesbrook Mall, Vancouver, B.C. Canada V6T 2A3 
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I. INTRODUCTION 


The determination of the phase structure of strongly interacting matter is one of the 
major challenges in theory and experiment. In our contribution we study the phase transition 
from hadrons to a quark-gluon plasma in asymmetric matter. Based on a two-phase model 
we describe new qualitative features that arise as a function of the isospin and that may be 
relevant for the phenomenological description of highly energetic heavy-ion collisions. 

The creation of the quark-gluon plasma and the elucidation of its properties is the fun¬ 
damental goal in relativistic heavy-ion physics. Although experiments at the Brookhaven 
AGS and at the CERN SPS have provided evidence that in highly energetic collisions an 
intermediate state is formed which is characterized by very high energy densities [[I], it is 
not clear if this new phase of matter has actually been reached (for a recent review, see @). 

On the theoretical side the determination of the equation of state (EOS) of strongly 
interacting matter at high temperatures and densities is a major problem in QCD. One of 
the most important issues here is revealing the phase diagram of matter and the nature of 
the transition. Most reliable results come from lattice calculations for systems at zero baryon 
density (for an updated status report see [|]). These calculations clearly indicate a first- 
order transition for pure SU(3) gauge theories 0. When dynamical quarks are added the 
situation is less transparent. The nature of the transition not only depends on the number 
of flavors but also on their masses ||||, particularly on the mass of the strange quark. For 
example, three light flavors show a first-order transition that is absent for larger masses ||. 
Despite the need of corresponding information in applications to heavy-ion collisions, lattice 
calculations at finite baryon density presently remain at an early stage without a definite 
conclusion [Q. 

To the best of our knowledge, there are no lattice calculations available for asymmetric 
systems, i.e ., systems at finite baryon density and with nonzero isospin. However, it has been 
demonstrated in case of the liquid-gas phase transition of nuclear matter that asymmetric 
systems have qualitatively different properties than symmetric systems . In particular, 
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the transition in an asymmetric system is continuous and of second order according to 
Ehrenfest’s definition, in contrast to the discontinuous, first-order transition in symmetric 
systems. The goal of the present work is to extend these studies to the deconfinement phase 
transition of asymmetric matter. Since lattice calculations cannot provide the required 
information at present, we resort to a separate description of the hadronic and the quark- 
gluon phase. This approach is not new and was previously employed to investigate the 
existence of quark cores in neutron stars |U|[TI|. Although such a formulation oversimplifies 
the full complexity of the problem, it is useful for providing a first orientation to the new 
features that arise when a system with two conserved charges, i.e., baryon number and 
isospin, undergoes a transition. 

The EOS for the quark-gluon phase is based on a simple “bag model”, where the thermo¬ 
dynamic potential is modeled as a sum of a perturbative contribution and a nonperturbative 
part which is represented by the bag constant. Along these lines we adopt a model which 
is based on two massless flavors and which includes the first-order correction in the strong 
coupling constant a s B0- 

The hadronic EOS is generated in a relativistic mean-field approach involving the inter¬ 
action of Dirac nucleons with isoscalar scalar and vector meson fields and with isovector p 


meson and pion fields [|13|,[14|. The central object in this approach is an energy functional 
of the meson mean fields, which, in principle, can be formulated without reference to an ex¬ 
plicit Lagrangian [[Uj. It has been shown in numerous applications that this type of model 
provides an excellent description of nuclear matter and of bulk nuclear properties through¬ 
out the periodic table jp~4| |22||. Thus we can calibrate our hadronic equation of state at 
zero temperature and normal nuclear densities, and then extrapolate into the regime of high 
density and temperature appropriate for the phase transition. The model we study involves 
cubic and quartic isoscalar self-couplings. Although nonlinear p meson interactions are pos¬ 
sible |23[], the nuclear symmetry energy is achieved using the simplest possible coupling of 
the p meson held to the nucleon. The pions are incorporated by adding their one-body 
contribution, i.e., the contribution of an ideal Bose gas with an effective chemical potential, 
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to the thermodynamic potential. The effective chemical potential is chosen such that the p 
meson held couples to the total isospin density, which receives a contribution from nucleons 
and pions. 

At low temperatures, the qualitative features studied here depend strongly on the nuclear 
symmetry energy. At higher temperatures, the inclusion of the pions is crucial due to the 
onset of a pion condensate in very asymmetric systems. The pions are incorporated in a 
very simple manner. Generally the isovector dependence of the hadronic EOS may be more 
complicated than in our model; however, these properties are not well known, particularly 
when one is far from symmetric matter. In this context it is important to realize that the 
hadronic model is formulated in terms of effective degrees of freedom. The two conserved 
charges, baryon number and total isospin, are the relevant quantities in the thermodynamic 
analysis whereas the particle numbers of the individual species may change when parameters, 
e.g., temperature and pressure, vary. In the present case the pion condensed phase resembles 
closely matter containing A-resonances, which we could include in our model as well as other 
unstable nucleon resonances. However, we expect that this will not change the qualitative 
features of the EOS, since it depends only on the total baryon number and isospin. 

Although our EOS can be calibrated very accurately at zero temperature and normal 
nuclear densities, the extrapolation into the regime of high densities is not unique. In the 


context of neutron star calculations it has recently been demonstrated [g3j that the high 
density EOS is sensitive to combinations of parameters that are difficult to calibrate at nor¬ 
mal densities. As a consequence, large uncertainties in the predictions arise. The density 
regime where the transition of hadronic matter to a quark-gluon plasma takes place corre¬ 
sponds to the interior region of a neutron star. Therefore similar difficulties can be expected. 
Although this limits the quantitative results of our analysis, the qualitative features of the 
basic physics remain unchanged. 

Studies of the deconfinement phase transition in terms of a separate description of the 
two phases have a long history in the literature. In particular, in the context of neutron 


star matter, many calculations similar to ours have been performed earlier |q.|Td| [i3| . |24|.[2o 
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The first analysis which considered the transition in a system with two conserved charges 
was performed by Glendenning [Q, with the restriction to zero temperature and to neutron 
star matter, i.e., matter in beta-decay equilibrium. Along these lines the effect of neutrino 
trapping on the composition of cold protoneutron stars with quark cores has been investi¬ 
gated more recently [J23J. We will extend the analysis to systems at arbitrary temperature, 
baryon density and isospin. By definition, the two-phase model cannot reveal the actual 
nature of the transition. In contrast, other phenomenological approaches which are directly 
inspired by the symmetries of QCD, for example linear sigma models |26|, have been applied 
to shed more light on the fundamental properties of the transition. However, whether it is 
of first or second order or a cross-over phenomenon can only be established by full QCD 
studies, i.e., by lattice calculations. We believe that, until this information is available, our 
approach is useful for providing a first orientation of the qualitative features which arise 
in the phase diagram of strongly interacting matter at finite baryon density and isospin. 
At zero baryon density the two-phase model implies a first-order transition. Starting from 
there our results are qualitatively similar to what was found in case of the liquid-gas phase 
transition of nuclear matter [j.9|]. Primarily, the transition of hadrons to quarks and gluons 
in asymmetric systems is of second order, in contrast to the first-order transition that oc¬ 
curs in symmetric matter (which behaves like a one-component system). That means that 
some thermodynamic variables, e.g., the entropy, are continuous and smooth throughout 
the transition, whereas the first-order transition is characterized by discontinuities in these 
quantities. Furthermore, we find that the transition regime in terms of densities and tem¬ 
peratures varies significantly with the isospin of the system. In particular the density at 
the onset of the transition is always lower in asymmetric matter. This result led us to the 
conclusion that it is more likely that the transition region is reached in collisions of very 
neutron-rich nuclei. 

The outline of this paper is as follows: In Sec. II, we briefly summarize the general 
discussion of phase transitions in multicomponent systems which can be found in greater 
detail in Ref. ||. We focus on the new aspects which arise when a thermodynamic sys- 
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tem is described by two different equations of state. In Sec. Ill, we describe the models 
for the two phases and summarize the relations that determine the equation of state. In 
Sec. IV, we apply our model to the deconfinement phase transition. We study the phase 
coexistence region and illustrate different thermodynamic processes by the aid of various 
phase diagrams and Maxwell constructions. We consider special situations which arise dur¬ 
ing ultra-relativistic heavy-ion collisions. We also discuss the parameter dependence of the 
thermodynamic features. Sec. V contains a short summary. 

II. PHASE TRANSITIONS IN MULTICOMPONENT SYSTEMS 

We consider a system characterized by a set of n mutually commuting charges Qi. The 
thermodynamic features of such a system have been discussed earlier [:7j.|9[] in connection 
with the liquid-gas phase transition of asymmetric nuclear matter. Here we will focus on 
the situation where the system can exist in two phases each characterized by a different 
equation of state. The partition function might be generated by a Hamiltonian which is 
formulated in terms of different degrees of freedom in each phase [|8|. For example, the 
elementary constituents (quarks and gluons) would be used in the first phase and composite 
particles (hadrons) in the second. As mentioned in the introduction, it is important to keep in 
mind that the conserved charges are the relevant quantities in the thermodynamic analysis. 
A conserved charge does not necessarily imply an independent particle species which may 
change during a process (by decay, for example); it includes any conserved quantity resulting 
from the underlying symmetries of the system, for example, isospin, baryon number, electric 
charge, etc. 

To describe the equilibrium state of the system in phase a (a = 1,2) enclosed in a volume 
V we choose the Helmholtz free energy F 

F a (T, V, Qi) = VF a (T, pi) , a = 1,2 (1) 


with 
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Qi_ 

V 


i — 1,... ,n . 


Pi = 


Depending on the value of the free energy the system will be realized in one phase or the 
other or in a mixture of both, indicating a phase transition. To be more precise, the system 
in phase a will be stable against separation into two phases if the free energy of the single 
phase a is lower than the free energy in all two-phase configurations. This requirement can 
be formulated as 

Fot(T, pi) < (1 — A)^ 7 a (T, pf) + \Tg(T, pf) , a,/? =1,2, (2) 

with 

Pi = (1 - A )p? + Apf , 0 < A < 1 . (3) 

Note that this also implies 


Fa(T,pi) < Tp{T,pi) , 


( 4 ) 


so if Eq. (H) is satisfied, phase a is not only stable against phase separation, it is also the 
energetically favorable single-phase configuration. We assume that each single phase by 
itself describes a stable configuration^, so that Eq. (|2|) is satisfied for a — (3, he., the free 
energy in each phase is a convex function of the densities [^7| . The last equation (J3j) ensures 
that the overall charges are conserved: 


V Pi = V a p a { + with V = V a + . 


( 5 ) 


Equation @ is a global criterion for the stability of a one-phase configuration. Whenever 
it is violated, a system with two phases is energetically favorable. The phase coexistence 


1 When it actually happens that the free energy of one of the single phases describes an unstable 
configuration, e.g., due to the liquid-gas phase transition in nuclear matter, we take T a to be the 
free energy resulting from a Maxwell construction in the unstable region, which is convex. 
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is governed by the Gibbs’ conditions demanding equal pressure p, temperature T and 
chemical potentials /i, in both phases 



( 6 ) 


(7) 


It is important to realize, however, that under conditions of phase separation, Eq. (|2|) may 
still be valid locally , that means, in some small region U = {T,pi} of parameter space, but 


it may nevertheless be possible to find significantly different densities pf and pf in a larger 


domain U C G that violate this condition. This leads to the existence of metastable states. 


The two sets of densities {p“,pf} that satisfy Eqs. (|6|) and ( 0 ) form a surface in the 


parameter space {T, /?*}; this is the phase separation boundary, or binodal. For n conserved 
charges and two coexisting phases, Gibbs’ phase rule implies that the binodal is an n- 


dimensional surface It can also be shown that this surface encloses all points that lead 


to a single (unstable) configuration with a higher value for the free energy @ . 


The binodal surface determines the stability boundaries of the system. The mixed phase 
inside must be determined by a Maxwell construction. This is achieved by solving Eq. ([|) 
|7]|| for given values of p*, with pf and pf lying on the binodal surface. The free energy in 
the transition region is then given by 


HT,Pi) = (1 - A )F a {T,p?) + AJEg(T, pf) . 


( 8 ) 


Densities related to other extensive quantities can be computed accordingly. 

These ideas are illustrated in Fig. ([!]). We show the free energy density in both phases 
as a function of the density coordinates {p*} at constant temperature. At the location of 
the points B i and B 2l the system enters the binodal region. In between, a phase mixture is 
energetically favorable. The free energy of the mixed phase as it results from the Maxwell 
construction Eq. (§) is indicated by the solid line which connects B\ with B 2 . The dotted 
lines correspond to the free energy of the single-phase configurations. Between the point of 
intersection / and B\ and between / and B 2 lie the metastable states of phase 1 and phase 2, 
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respectively. The configurations beyond / are unstable. Here not only a phase mixture but 
also the single (metastable) configuration of the opposite phase is energetically favorable. 
Instability normally arises due to fluctuations, which no longer can be restored by the system. 
This is signaled by the violation of a local stability criterion, e.g., a negative compressibility. 
However, our description is based purely on the energetics and cannot provide information 
about the physical mechanism which leads to instability. This is certainly an indication that 
the separate description of the two phases breaks down at this point. 

We close this section by specializing the general formalism to asymmetric matter. In the 
hadronic phase it consists of strongly interacting nucleons and mesons, and in the decon- 
fined phase, weakly interacting quarks and gluons. For two quark flavors such a system is 
characterized by two conserved charges: the total baryon number 

N b = V ps (9) 


and the third component of isospin 


h= ^Vps ■ 


( 10 ) 


Thus we have for the thermodynamic potential, or equivalently the pressure, 
~^( T > M%> Ms) = ~P a = F a { T \ PB , Ps) - MbPB - IjMsPs , 


( 11 ) 


where the baryon and isospin chemical potentials are given by 

' dTg \ 

. ®Pb ) t , P3 


Mb = 


, A*3=2 


0Fa \ 

/ T,p B 


( 12 ) 


In the discussion of asymmetric systems it is also useful to introduce the isospin ratio 

h Ps 


x = 


Nb %Pb 

so that the free energy can be rewritten as 


(13) 


Fa(.T,p B ,p 3 ) = F a (T,p B ,x) . 


(14) 
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III. THE EQUATION OF STATE 


To generate the nuclear equation of state in the hadronic phase, we adopt the relativistic 
mean-held approach of Refs. |15|.23j involving valence Dirac nucleons and effective mesonic 
degrees of freedom, which are taken to be neutral scalar (0) and vector fields (V M ), plus the 
isovector p meson (b M ) and pion held (7r). The basic quantity here is the thermodynamic 
potential as a functional of the meson mean helds which, in principle, can be formulated 
without reference to an explicit lagrangian. For homogeneous nuclear matter, at finite 
temperature and chemical potentials (/i B ,/i 3 ), it is of the general form 


y(T, (jib, /U! 0, K, K, ir) = ^ m 2 (f> 2 - i ^ m 2 b M -b^ + ^ rn^ir 2 


AV(T, /x B , /i 3 ; 0, V0, b„, n) . 


(15) 


The hrst two terms in Eq. (|I5D, fl N /V and Vl M /V , are the one-body contributions due to va¬ 
lence nucleons and mesons. The nonlinear potential AV represents the unknown part of the 
thermodynamic potential, which includes the effects of nucleon exchange and correlations, 
as well as contributions from the quantum vacuum JI5|J3^ . 


Following Ref. [|J the one-body fermionic contribution can be written as 
Qjv 1 


V 


37T 2 - 


H 5 (v p ,M*) + H 5 (v n ,M*) 


(16) 


where the integral H 5 as well as others will be defined shortly. The corresponding contribu¬ 
tions of the heavy meson helds (0, U M , b M ) to 0^ are negligible in the relevant temperature 
range (T < 200 MeV). Thus we keep only the pions and write 


ft 


M 


V 


ftyr 

V 


6n 2 - 


B 5 (v n ,m n ) + -. 85 ( 0 , 771 *.) 


(17) 


where the hrst term is the contribution of the charged pions |p0| . |31|| , and the latter arises 
from the 7T°. The nonlinear potential AV can be expanded in a Taylor series in terms of 
the meson mean helds ||T5|,^3j. In the normal phase of nuclear matter, he., without a pion 
condensate, the pion mean held vanishes and the pions contribute only to the coefficients of 
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this series via loops. In practice, the series must be truncated, and the unknown coefficients 
serve as model parameters, which are chosen to reproduce certain empirical properties of 
equilibrium nuclear matter, as discussed below. For our calculation we choose the explicit 
form 

AV = - |^(VW") 2 . (18) 


which includes a subset of the meson self-interactions up to fourth order in the fields. In 
principle, terms which couple different fields and also nonlinear interactions involving the p 
meson held are allowed |15|. We also disregard contributions involving the pion mean held 
which arise in the pion condensed phase. Although disregarding these couplings is “unnat¬ 
ural”, the present model defined by Eq. (|l^) is already general enough for our purposes. 
First of all, it can be accurately calibrated at normal densities, and it can be related to the 
most common models discussed in the literature. Unless one assumes scalar self-interactions 
only, the nonlinear couplings cannot be uniquely constrained by the calibration procedure, 
leading to uncertainties in the high-density EOS. In particular, little information is available 
to constrain contributions which explicitly depend on the pion held. Moreover, among these 
nonlinear terms, the quartic vector coupling in Eq. (|T^) is the most important one at high 
densities ^3[. We will use this coupling to examine the uncertainties which arise in the 
quantitative predictions for the transition from hadrons to a quark-gluon plasma. 

The mean helds are determined by extremization of the thermodynamic potential 
Eq. (p~5D . The specihc form of AV, introduced in Eq. (|T8|) , together with Eqs. (pT|) and 
(□>. leads to the self-consistency equations 


9 2 s 2 gl 6 gt 


& = Ps, 


W 1 + 


9l C 

ml 6 


W 2 ) = 


9v 


m 


P > 


r> _ 9p 

2yU 


(19) 

( 20 ) 
( 21 ) 


where the scalar and baryon densities are given by 
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M* 

Ps — 2 

7T Z 

1 r 

PB ~ — 
n z L 


H 3 (u p ,M*) + H 3 (u n ,M*) 
t?3( u p , M*) + M*) 


The isospin density receives contributions from nucleons and pions 


( 22 ) 

(23) 


\pi = \p” + \pI (24) 

= -^2 (u p , M*) — C'nOv f t", 1 ] + nijr) . 

The scaled meson fields are $ = g s 0, PE = g v Po, and R = g p b 0 , with b 0 the timelike, neutral 
part of the p meson field. 

The baryon effective mass and effective chemical potentials are defined in terms of the 
meson mean fields and the baryon and isospin chemical potential as 


M* = M - d> , 

Up = p B + y - W ~ i R , 
v n = PR - y - w + ^R. 


The effective chemical potential for the pions is given by 


Uff — Up u n ps R , 

so that the p meson field in Eq. (^l|) couples to the total isospin density in Eq. 
also introduce the required integrals over the thermal distribution functions as 

POO / \ \ \ 

G n (p, M ) = J o k dk (j +e 0 [E(k,M)~lA - 1 + e P[E(k,M)+lA J ’ 


H n (p,M)= [ 
Jo 

for fermions and 


00 k n l dk 


+ 


o E(k,M)\ l + e Pl E ( k , M )-n] l + e p[E(k,M)+iA J ’ 


A n (p, m) = J k n 1 dk ^ 

B n (p,m) — y e p[E(k,m)-n] _ l 1 e /3[E(k,m)+n] _ l J ’ 


e f3[E(k,m)—fi] ^ g p[E(k,m)+n] \J ’ 

°° k n -'dk t 1 1 


(25) 

(26) 
(27) 


(28) 

We 


(29) 

(30) 


(31) 

(32) 
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for bosons, where E(k, M ) = ( k 2 + M 2 ) 1 / 2 , and n > 0 to ensure convergence is understood. 
Moreover, the boson integrals are subject to the constraint |//| < m. 

By inserting the solution of Eqs. ([T9|) —(pTD into Eq. ([15]), it is straightforward to compute 
the pressure 


^ frp T\J*\ frp 

P = -77 = —— [T, i/p, i/ n , M ) - — (T, z/ w 


V 


V 


V 


+ w2 + 4 14/4 + T2 r2 - wh $2 - TTs $3 - 

29? 24 2 gi 2gi 69 ? 24<rf 


(33) 


Other quantities like the free energy, the entropy, etc., follow by using standard thermody¬ 
namic relations. 

For the quark-gluon phase we adopt a bag-model type EOS |T2| involving massless u and 
d quarks (Nf = 2) 


y (T, Pb, P 3 ) - // 3 ) + b , 


(34) 


where Q peit /V is the perturbative expansion of the thermodynamic potential 



and b is the difference between the energy density of the perturbative and the nonpertur- 
bative QCD vacuum, he., the bag constant. The chemical potentials of the two flavors are 
connected to the baryon and isospin chemical potentials by 

1 1 11 

Pu — ~PB + ~Pz , Pd — ~PB - ~P3 , (36) 


which allows the computation of the corresponding densities according to 


Pb — ~ 


' d 

^V) T , 3 


:P3 = - 


' 8 Q' 


(37) 


T,hb 


We close this section with a short description of the calibration of the EOS. The model 
for the hadronic phase has six free parameters. Following Ref. [|15| we choose the value of 


14 
















the coupling £ and determine the other parameters so that the five equilibrium properties 
of nuclear matter, as listed in Table ill, are reproduced. The free parameter ( is chosen 
within the natural range 0 < ( < 0.06 To fix the two parameters in the quark- 

gluon phase, namely the bag constant b and the strong coupling constant a s , we proceed 
as follows. We specify the bag constant to b = 160 MeV/frn 3 and determine a s so that 
our model reproduces a transition temperature of 150 MeV at zero baryon density, which 


is in the currently accepted temperature range |34|]. As a consequence of this calibration a s 
depends on (. This dependence, however, is extremely weak. The value a s = 0.349 obtained 
for ( = 0 changes only by 0.1% when ( is varied within the natural range. 


IV. PHASE STRUCTURE OF ASYMMETRIC MATTER 

In this section we apply the formalism of Sec. II to the deconfinement phase transition. 
Within our simple two-phase model we consider highly excited asymmetric matter created 
in an ultra-relativistic collision of two heavy-ions. The dynamical evolution of the system 
during the collision is a very complex process and it is convenient to divide the reaction 
into three stages: the formation of a highly excited matter system, its expansion and the 
eventual decay. The question of whether the required energy densities can be achieved in 
an experiment to form a quark-gluon plasma during the first violent stage is a difficult one, 
which we will not attempt to answer here. We simply assume that such a state arises and 
ultimately reaches equilibrium at some finite temperature, density and pressure. We follow 
the subsequent expansion of the system during the second stage and concentrate on the new 
features of the phase diagram that arise as a function of the total isospin. Although our 
general discussion will encompass the entire range of x < 0, we will study systems with 
—0.2 < x < 0 and x = —0.5 in more detail to obtain estimates for the size of the new 
effects. The former range covers the isospin ratios which are experimentally accessible in 
heavy-ion collisions, e.g., £au+Au, aqj+u ~ —0.1, whereas the latter value corresponds to 
neutron matter relevant in astrophysical applications. 
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As discussed in the context of the liquid-gas phase transition of asymmetric nuclear 
matter |7].[3[], the phase structure of binary systems is more complex than in one-component 
systems. The basic quantity in the analysis is the phase separation boundary or binodal, 
which forms a two-dimensional surface in parameter space as discussed in Sec. II. Before we 
discuss the binodal structure of the deconfinement phase transition as it arises in our model, 
let us briefly point out the role of the pions. To this end, assume our system is prepared in 
the hadron phase at some fixed isospin ratio x, as defined in Eq. fll3D , and let us follow an 
isothermal compression. At some critical baryon density p c B the chemical potential of the 
pions will approach the (negative) value of the pion mass, i.e., 

v* - m *■ for p B -> p c B , (38) 


indicating the onset of a pion condensate |29|. The critical density depends on x and on T, 
and will be smaller in more asymmetric systems. At this point a macroscopically occupied 
7 t~ mode arises which is characterized by a nonvanishing pion mean field. In principle, 
the corresponding contributions to the thermodynamic potential in Eq. (0 must now be 
kept explicitly, and the value of the mean field is determined by minimization [RD1RS|. As 
discussed in the last section we will disregard the explicit mean-held contributions of the 
pions in the following. That means we neglect the pion-nucleon as well as the pion self¬ 
interactions. To incorporate the condensed phase consistently we will be satisfied to add the 
contribution of the zero-momentum states of the pions to the isospin density in Eq. (p 4|) 


\pl = —2m n K 2 + ■^—A 3 (—m n , m,) , 

Z Z7T z 


(39) 


where the effective chemical potential 


Vn = v p -Vn = 


(40) 


is held fixed in the condensed phase. The parameter n 2 is the continuous order parameter of 
this transition, which specifies the amount of isospin charge carried by the zero-momentum 
states fl30ljni. In contrast, the condensate does not contribute to the pressure of the system; 
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thus, one simply substitutes the value of u n , as given by Eq. fl40|) , in the expression for the 
hadronic pressure Eq.([33|). We will leave the study of the contributions arising from pionic 
interactions as an important topic for future work. 

We now return to the general discussion of the phase diagram in our system. We begin 
the analysis with ( = 0 for the value of the quartic vector meson coupling in the hadronic 
EOS. We will study the influence of this parameter at the end of this section. 

The binodal is determined by Gibbs’ conditions (|J) and (0) which enforce equal pressure 
and chemical potentials for the hadronic phase and for the quark-gluon phase in equilibrium. 
In general these two phases have a different isospin ratio x. The binodal surface in {x,p,T} 
space is indicated in Fig. Q. A sequence of several slices at fixed T is shown. For a given 
temperature, the binodal section is divided into two branches. One branch, at the lower 
pressure, describes the system in the hadron phase, while the other branch, at the higher 
pressure, describes the quark-gluon phase. These two branches contain the beginning and 
ending configurations of the phase transition. The shape of the binodal slices changes 
drastically between T = 0 and T = 150 MeV. One observes that the enclosed area becomes 
smaller with increasing temperature until, at T c = 150 MeV, the two branches coalesce 
to a single line at constant pressure p = 19.08 MeV/fm 3 with isospin ratios ranging from 
x ~ —4.5 to x — 0. This line of critical points marks the transition point of systems at zero 
baryon density. 

Another important observation is that the pressure steadily decreases when the tem¬ 
perature increases. At the same time, the transition regime, i.e., the difference in pressure 
between the two phases at the beginning and at the end of a transition, becomes smaller. 
Also indicated in Fig. (||) are the special configurations that separate into two phases with 
the same isospin ratio. The coordinates of these configurations form the two lines of equal 
concentration (LEG), one of which coincides with the projection x — 0 in our system. The 
isospin ratio of the second LEG increases with temperature and eventually the two LEC 
intersect at T « 137.2 MeV. For a fixed temperature the pressure necessarily attains an 
extremum at the points of equal concentration B 
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Generally, the hadronic configurations are divided into normal states and states exhibit¬ 
ing a pion condensate. The coordinates of the configurations at the onset of the condensation 
also form a two-dimensional surface. This critical surface can be parametrized in terms of 
two coordinates in {x,p, T} space, e.g., as p c (x,T). The intersection of the critical surface 
with the binodal forms the line of three phases (LTP), which marks the onset of the pion 
condensate on the hadronic branch^. In Fig. (||) these states are emphasized by the bold 
dashed lines. The isospin ratio of the LTP decreases with increasing temperature. Above 
T fa 136.2 MeV the critical pressure curve and the binodal no longer intersect. 

The binodal in {pb, x,T} space is indicated in Fig. (|3]), where several branches at different 
x are projected onto the {pb,T} plane. The onset of the condensate, i.e., the critical 
temperature as a function of the baryon density and isospin ratio, is marked by an open 
circle. Note the change of density in the two-phase region. The most extreme situations are 
encountered at T = 0. For example, in symmetric matter the density increases between the 
onset and the completion of the transition by more than 3.5 p° B . This density regime becomes 
even larger in asymmetric systems. Another important observation is that the density at 
the onset decreases with increasing asymmetry. If we compare these points for x = 0 and 
x = —0.2 at T = 0 we find an decrease of A p B ~ 0.5p%. 

The role of the pions becomes more apparent in Fig. (}4]). Part (a) indicates one slice 
of Fig. (|H) at T = 50 MeV. Also included is the projection of the critical surface in the 
form of the critical pressure curve p c (x). Nuclear matter in the normal phase resides on the 
right-hand side, and the states with a pion condensate on the left-hand side. The critical 
curve intersects the binodal at TP. In contrast, the pion contributions are not included 
in part (b). With nucleons only, the isospin ratio in the hadronic branch is limited by 

2 Since there is no two-phase coexistence in a Bose condensation, the LTP is not a line of triple 
points which could arise in a binary system . The LTP is merely a boundary between areas in 
parameter space which describe two different two-phase configurations. 
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—0.5 < x < 0 which leads to a widely open binodal surface. Furthermore, only one point of 
equal concentration arises (at x — 0). 

To illustrate the phase-separation scenarios, we study the behavior of matter under an 
isothermal compression. Assume that the system is initially prepared in the hadronic phase 
with isospin ratio x = —0.4 and p < 50 MeV/frn 3 . The situation is also indicated in part 
(a) of Fig. (|. During the compression, the system first crosses the critical pressure curve at 
the point A, where pion condensation begins. Starting with k 2 = 0 at A the contribution of 
the condensate to the (fixed) isospin ratio increases during the compression. At the point 
B the two-phase region is encountered and now a quark-gluon phase is about to emerge. 
The solution of Gibbs’ conditions (||) and (0) determine the density and the isospin ratio x c 
of this new phase, which occurs at the point labeled C. As the system is compressed, the 
total isospin ratio x remains fixed, as dictated by the conservation laws, but the hadronic 
phase evolves from B to E , while the quark-gluon phase evolves from C to D. At the point 
D, the system leaves the region of instability. The original hadronic phase (which now has 
no pion condensate) is about to disappear, and it exists in an infinitesimal volume with a 
density and isospin ratio x E corresponding to the point E. 

An important observation is that if the system encounters the binodal between the two 
points of equal concentration then the quark-gluon phase is always more asymmetric and 
has a higher baryon density than the hadronic phase. This behavior is due to the symmetry 
energy, which is higher in the hadronic phase. Therefore it is energetically favorable for 
asymmetric matter to separate into a denser and more asymmetric quark-gluon phase and 
a hadronic phase that is more dilute and symmetric. 

The energetics also implies a significant simplification if the system is initially prepared 
with an isospin ratio corresponding to one of the points EC. It is now energetically favorable 
for the system to separate into two phases with the same composition. In this so-called 
“indifferent equilibrium” matter behaves like a one-component system (see also the 
discussion in 0). In contrast to the previous case the system becomes unstable at one of 
the points EC in Fig. [| and stays at this point until the transition is completed. Note that 
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the location of the point EC at x = 0 is model independent since symmetric matter is always 
the configuration with lowest energy. The location of the second point depends on specific 
model features, particularly on the predicted values for the symmetry energy and on the 
pionic contributions. 

The results of the corresponding Maxwell constructions are indicated in Fig. [5]. We 
consider different isospin ratios at T = 0 in part (a) and at T = 50 MeV in part (b), 
respectively. For symmetric matter (x — 0), we obtain the familiar result with a constant 
vapor pressure, represented by a horizontal line. At x < 0 the equal demand of Gibbs’ 
conditions and the conservation of the overall isospin forces the pressure to change during the 
transition. Also shown in Fig. [5| are isotherms for “neutron matter” (x = —0.5). The name 
is somewhat misleading since the system contains a finite number of protons even without 
the pions. ft is common in the literature 0HO-0 to treat neutron matter as a single¬ 
component system. The qualitatively different features which arise in a proper treatment 
as a system with two conserved charges have been pointed out earlier by Glendenning || 
in the context of neutron star calculations (see also Ref. ]23j]). In addition, part (b) of 
Fig. [5] also includes a system with an isospin ratio which corresponds to the second point of 
equal concentration. As discussed above, this system also stays at a constant vapor pressure 
throughout the transition. 

The change of the pressure throughout the phase separation in asymmetric systems is 
an indication of a smoother transition than in symmetric or one-component systems. A 
more precise characterization of the transition can be obtained from the isobaric process 
indicated in Fig. ([]). The entropy as a function of temperature for three different isospin 
ratios is shown. In symmetric matter we find a discontinuity at the transition temperature, 
which gives rise to a latent heat Ql, i.e., all the heat is used to convert hadrons into 
quarks. In asymmetric matter the entropy remains smooth but the temperature now changes 
throughout the transition so that no strict latent heat can be assigned to this process. For 
example, the isospin ratio x = —0.2 shown in Fig. (;6!) gives rise to a temperature change 
of AT fa 30 MeV. As discussed in more detail in Ref. |J, these are the features which 
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lead to a first-order transition in symmetric matter and to a second-order transition in the 
asymmetric system, according to Ehrenfest’s definition. 


The processes we have considered so far correspond to very specific idealized situations 
and served mainly to study the basic thermodynamic features of our system. Generally, it 
is not clear if the expansion of the highly excited matter created in a heavy-ion collision 
can be described by a distinct thermodynamic process. However, the situation simplifies at 
extremely high collision energies. In this so-called scaling regime, the colliding nuclei become 
transparent and the excited plasma is created between the two nuclei when they recede from 
the collision point. The basic assumption here is that the total entropy is created in the 
initial stage of the collision and remains nearly constant during the later evolution. Under 
this assumption, the expansion of the plasma is nearly adiabatic and can be described by 
relativistic hydrodynamics |]37| , |38| . 


To make contact with this discussion, we also studied adiabatic processes. The properties 
of the system in the (p B , T ) plane can be studied in Fig. (|7|). Adiabats for various values of 
the entropy per baryon are indicated for symmetric matter in part (a) and for asymmetric 
matter in part (b), respectively. The dashed lines are the coexistence curves, i.e., the 
projection of the hadronic and the quark-gluon branch of the binodal, which intersect at 
the critical point CP. The curves in the coexistence region are the result of a Maxwell 
construction. Note the behavior of the system in the mixed phase. Throughout the transition 
the temperature decreases with increasing density, in contrast to the pure phases where the 
temperature increases monotonically. This peculiar behavior arises because for a given 
baryon density and temperature, the entropy in the quark-gluon phase is always higher. 
If we think in terms of a “boiling” or “condensation” process and compare to a liquid-gas 
transition one would expect the opposite behavior []. In an ordinary liquid the high-density 


3 A system which exhibits a similar behavior is the liquid-solid phase transition of 3 He, for which 
the entropy in the high-density solid phase is higher than in the low-density liquid phase due to 
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phase has the lower entropy and the temperature would rise steadily during the transition. 

If we compare the extension of the transition regime in the symmetric and asymmetric 
systems, we find a similar quantitative difference as observed in the previous discussion. 
Generally, for a given value of the entropy per baryon, the onset of the transition occurs 
at a lower density and temperature in asymmetric matter. For example, if we compare the 
adiabat S/Nb = 2 in Fig. and (0b), the density decreases by A ps ~ p^/3 and the 
temperature by AT & 4 MeV. Compared to the overall scale set by the critical temperature 
of 150 MeV this temperature difference is certainly less significant. However, the qualitative 
difference between the two systems could be an indication that the transition region is 
experimentally better accessible in a collision of extremely neutron-rich systems, which might 
be obtainable with radioactive beams. 

Finally, we briefly discuss how the properties of the phase transition change when the 
parameter ( in the hadronic EOS is varied. This point has been discussed in Ref. p3[ for pure 
neutron matter and here we will extend the analysis to arbitrary isospin ratios. Up to now 
the hadronic EOS in Eq. ( |55| ) was generated without the quartic vector meson interaction. 
In the following we allow a nonzero coupling within the natural range 0 < ( < 0.06. 

Figure (^j) shows the baryon density at the onset and at the end of the transition for 
the isospin ratios x = 0, —0.2, —0.5 at T — 0. Note that all values of ( reproduce identical 
properties of nuclear matter and a transition temperature of T — 150 MeV at zero baryon 
density. Generally, one observes that for arbitrary isospin ratios, increasing ( softens the 
EOS, which leads to substantially higher transition densities. The increase in the symmet¬ 
ric system is remarkable (note the logarithmic scale); at sufficiently large values of ( the 
transition vanishes altogether. The transition is driven by the energetics in both phases and 
therefore this behavior is most easily understood by examining the high-density limit of the 
energy per baryon |23|]. Using Eq. (|35|) and the relation £ = 3p + 46/3, one obtains for the 


the spin disorder present in the solid phase [39] 
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quark phase 

,,}}-L £ q/ Pr = c Q p z 3 ' ( 41 ) 

The asymptotic behavior of the hadronic EOS depends on £ and x. For x < 0, ( > 0 or 
x = 0, £ = 0 the quadratic terms in Eq. ([Tf) dominate, which leads to 


lim £ h /pb oc p B ■ (42) 

pB^OO 

Thus, at sufficiently high densities, the hadronic matter always has higher energy compared 
to the quark phase and a transition is possible. 

A special situation arises for x = 0 and £ > 0. In this case the leading behavior is, up to 
the prefactor, identical to Eq. (]4l|) , he., 


lim Eh/pb = chPb 3 • 

PB ^C50 


(43) 


By comparing the factors Cq and Ch in Eq. (|4]) and in Eq. (|43|) , respectively, one arrives at 
the remarable result that a transition is possible only if (see also the discussion in 


4 \ 1/3 


> + (^) 


> 1 + 


2 a s 

37T 


(44) 


This implies that for sufficiently large values of the nonlinear coupling, symmetric matter 
remains in the hadron phase. Note that we would arrive at this result for arbitrary isospin 


ratios if we included an additional quartic p meson coupling. In accordance with Ref. |23 


we conclude that the parameter dependence of the hadronic EOS at high densities leads to 
large uncertainties in the predictions for the transition region of the deconfinement phase 
transition (if it exists). In principle, reproducing the quantitative features of the transition 
could be a tool for calibrating the hadronic EOS better in the high-density region. However, 
more information is required than available at present, and quantitative predictions are 
limited. 
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V. SUMMARY 


In this paper we studied the deconfinement phase transition from hadronic matter to 
a quark-gluon plasma as a function of the isospin. The analysis was based on a separate 
description of the hadronic and the quark-gluon phase. We have demonstrated that in a 
consistent treatment of the two conserved charges, baryon number and isospin, qualitatively 
new features arise. Most importantly, the phase transition in asymmetric matter is generally 
continuous and of second-order. This is in contrast to the discontinuous first-order behavior 
of symmetric systems, which arise as singular points in the phase diagram of asymmetric 
matter. 

The analysis was based on the determination of the phase separation boundary, the bin- 
odal, which indicates the region where a separation into two phases is energetically favorable. 
The binodal is determined by Gibbs’ criteria for phase equilibrium. In a system with two 
conserved charges, it forms a two-dimensional surface in parameter space, in contrast to the 
one-dimensional surface in one-component systems. The greater dimensionality leads to the 
qualitatively new behavior. Primarily, the transition in asymmetric matter is continuous 
and of second order. We saw that the pressure, temperature and isospin ratios in the partic¬ 
ipating phases change throughout the transition. That means that in general, asymmetric 
systems pass through a transition region rather than staying at a transition point which 
leads to discontinuities in one-component and symmetric systems. 

To apply these results to the deconfinement phase transition in nuclear matter, we em¬ 
ployed a simple two-phase model. The hadronic EOS was generated in a relativistic mean- 
field model involving the interaction of baryons with isoscalar scalar and vector fields and 
with the isovector p meson and pion field. This model involves cubic and quartic isoscalar 
self-couplings and allows for an accurate calibration at normal nuclear densities to reproduce 
bulk properties of finite nuclei and nuclear matter. For the description of the quark-gluon 
phase, we adopted a bag-model type EOS involving massless u and d quarks. In this frame¬ 
work thermodynamic quantities consists of two parts. A perturbative contribution, where 
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we included the first order correction in the strong coupling constant, and a nonperturbative 
contribution, represented by the bag constant, which accounts for the QCD trace anomaly. 
Guided by recent results from lattice calculations, we chose the bag constant and the strong 
coupling constant to reproduce a transition temperature of 150 MeV for matter at zero 
baryon density. However, the calibration of the hadronic EOS at normal densities is not 
unique and predictions at high densities are sensitive to the model parameters. By studying 
variations of the quart ic vector meson coupling within the natural range, we found signifi¬ 
cant uncertainties in the predictions for the onset of the deconfinement phase transitions in 
systems at finite baryon density. For sufficiently large values of this coupling the hadronic 
EOS becomes so soft that the transition vanishes altogether. 

This two-phase model was then used to study the deconfinement phase transition in 
highly excited matter that is produced in ultra-relativistic heavy-ion collisions. By con¬ 
struction, an approach based on a separate description of the two phases cannot reveal the 
actual nature of the transition. However, until more concrete and reliable lattice calculations 
are available, we believe that despite this limitation our analysis is useful for providing a 
first orientation and concrete description of the qualitative features which arise in the phase 
diagram of strongly interacting asymmetric matter. There are several significant differences 
between the phase diagram for an asymmetric system and that for symmetric matter. In the 
first place, the location of the transition region in parameter space depends on the isospin 
ratio of the system. Generally, the onset of phase separation occurs at lower baryon densities 
and temperatures in more asymmetric systems. More importantly, the dimensionality of the 
phase-separation region is larger in asymmetric matter leading to a continuous transition. 
This implies that the thermodynamic properties of the participating phases change through¬ 
out the transition. Although these effects are small for realistic isospin ratios x ~ —0.1, the 
trend in the isospin dependence to lower densities and temperatures could be an indication 
that the transition region is easier to reach in collisions of very neutron rich nuclei, which 
might be created in radioactive beam facilities. 

The analysis also revealed interesting physics in the isovector channel. Most importantly, 
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the incorporation of the pions into the model for the hadronic phase was crucial for the shape 
of the binodal surface and therefore also for the path of the system through the mixed-phase 
region. The salient feature is the onset of a pion condensate at large baryon densities and in 
very asymmetric systems. This phenomenon is certainly of minor practical importance since 
the relevant isospin ratios lie beyond the experimentally accessible region. However, inter¬ 
est in Bose condensation has recently been revived in the literature, particularly since the 
possibility of kaon condensation in neutron star matter has been proposed [[h|]. Analyzing 
these effects in connection with the deconfinement phase transition provides an important 
topic for future work on this problem. 
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TABLES 


TABLE I. Equilibrium Properties of Nuclear Matter 
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FIGURE CAPTIONS 


FIG. 1. The Maxwell construction for a system with more than one conserved charge. The free 
energy density as a function of the densities {p,} at constant temperature is shown. The Maxwell 
construction is indicated by the segment B\B- 2 . The course of the free energy in each single phase 
is denoted by the dotted lines. The metastable configurations are located on B\I and B-jI between 
the binodal and the point where the free energy curves intersect. 

FIG. 2. The binodal surface indicating the two-dimensional phase-coexistence boundary in 
(p, T, x ) space. A sequence of binodal sections at fixed temperature is shown. The points of equal 
concentrations (EC) and the points of three phases (TP) are indicated. 

FIG. 3. Binodal section in ( ps,T,x ) space. The projection of several binodal branches at dif¬ 
ferent x onto the ( pb,T) plane is shown. The isospin ratios are x = 0, —0.1, —0.2, —0.3, —0.4, —0.5, 
starting from the right. The open circles mark the critical temperature for the onset of the pion 
condensate. 

FIG. 4. (a) Binodal section at T = 50MeV including the contributions of the pions in the 

hadronic phase. The points B through E denote phases participating in a phase transition. At 
A the system crosses the critical pressure curve which marks the onset of pion condensation. The 
point of three phases (TP) and the points of equal concentration (EC) are also indicated, (b) 
Binodal section at T = 50MeV without the pionic contributions. 

FIG. 5. (a) Isotherms for different values of x at T = 0. The isospin ratios are 

x = 0,—0.1,—0.2,—0.3,—0.4,—0.5 from the top to the bottom, (b) Isotherms for different val¬ 
ues of x at T = 50MeV. The isospin ratios are the same as in part (a). In addition the isotherm 
with x EC = —1.44, which corresponds to the second point of equal concentration, is also indicated. 

FIG. 6. Entropy per baryon as a function of temperature at constant pressure for various 
isospin ratios. 
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FIG. 7. Properties of nuclear matter as functions of temperature and density. The dashed 
lines denote the coexistence curves (CE), which intersect at the critical point (CP). The solid lines 
indicate adiabats. The lines between the coexistence curves are the stable mixed-phase configura¬ 
tions. Part (a) shows results for symmetric matter (x = 0) and part (b), for asymmetric matter 
with x = —0.2. 

FIG. 8. Baryon density at the onset and at the end of the transition as a function of the quartic 
vector meson coupling (. 
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